Gravitational waves and pulsating stars: What can we learn from future observations? 
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We present new results for pulsating stars in general relativity. First we show that the so-called gravitational-wave modes of 
a neutron star can be excited when a gravitational wave impinges on the star. Numerical simulations suggest that the modes 
may be astrophysically relevant, and we discuss whether they will be observable with future gravitational-wave detectors. We 
also discuss how such observations could lead to estimates of both the radius and the mass of a neutron star, and thus put 
constraints on the nuclear equation of state. 
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Pulsating stars have been studied within Einstein's 
general theory of relativity for thirty years. Hence it 
is somewhat surprising that recent results have altered 
our understanding of these problems considerably. We 
now know that there are oscillation modes that are di- 
rectly associated with the curvature of spacetime. But 
it is not yet clear whether these modes are of any rele- 
vance to astrophysics, or if they are simply a curiosity of 
the mathematical theory. In order to demonstrate their 
relevance one must show that the modes are excited in 
a realistic astrophysical situation, such as gravitational 
collapse to form a neutron star. Unfortunately, we do 
not yet have the ability to perform such calculations in a 
relativistic framework. While waiting for full 3D numer- 
ical relativity to provide a definitive statement regarding 
the relevance of the gravitational-wave modes of a com- 
pact star, we must settle for suggestive calculations. The 
purpose of this letter is to describe one such calculation: 
We consider the scattering of gravitational-wave packets 
off a uniform density star. In this case we find that the 
gravitational-wave modes are indeed excited. 

Recent work on the pulsation properties of compact 
stars has shown that general relativity plays an impor- 
tant role: A relativistic star exhibits two distinct sets of 
pulsation modes. One is slowly damped and corresponds 
to the well-known fluid modes , while the second is 

rapidly damped and has no analogue in Newtonian the- 
ory. The new modes have been termed w-modes because 
they are closely associated with gravitational waves [Q. 
We have recently managed to put the w-modes in con- 
text and provide a better understanding of the physics 
involved |,§7|]. 

The equations describing perturbed nonrotating stars 
in general relativity split into two classes Q|. Polar 
perturbations correspond to zonal compressions of the 
star, whereas axial ones induce differential rotation in 
the fluid. In general, the polar problem corresponds to 
two coupled wave equations: one that represents the fluid 



motion and one for the gravitational waves Q . The fluid 
pulsation modes (eg. the f-mode) belong to the polar 
class of perturbations. Meanwhile, axial perturbations 
are described by a single, homogeneous wave equation 
||. The axial problem is thus considerably simpler. Nev- 
ertheless, it attained very little attention until recently. 
This is mainly because an axial perturbation can only 
induce pulsations in the stellar fluid through a non-zero 
shear modulus || , or through coupling to the polar per- 
turbations if the star is rotating jlO|. However, the axial 
problem is interesting also for the simplest stellar mod- 
els: If the star can be made ultracompact (R < 3M in 
geometrical units c = G = 1) the peak of the exterior 
curvature potential barrier (that is familiar from black- 
hole perturbation theory) will be unveiled. Then gravi- 
tational waves that impinge on the star can be trapped. 
That such trapped modes exist has been demonstrated 
|]lT| , ^2| , but they are unlikely to be of any great astrophys- 
ical relevance: A sufficiently compact star will probably 
never form. 

More importantly, one can argue that the axial prob- 
lem is relevant also for less compact stars (the canonical 
values R — 10 km and M — 1.4M Q for neutron stars lead 
to M/R « 0.2), especially as far as gravitational waves 
are concerned. For black holes, when the gravitational 
waves are the only active agent, the axial problem leads 
to a pulsation spectrum that is identical to the polar one. 
Knowing this, the results of a detailed mode-survey for 
uniform density stars do not come as a complete surprise: 
The axial and the polar w-mode spectra are remarkably 
similar Q. This is an important result because of its 
implications for physical interpretations of the w-modes. 
Since the axial modes do not couple to the stellar fluid 
one cannot invoke the fluid in an explanation of them. 
The w-modes must arise because gravitational waves can 
be temporarily trapped in the "bowl of spacetime curva- 
ture" provided by the mass of the star || . 

Are the w-modes excited?. — Our understanding of the 
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role of a dynamic spacetime for pulsating stars has been 
improved, but the most important question remains. As 
yet it has not been established that the w-modes can con- 
tribute to observable gravitational waves, and thus play a 
role in astrophysics. The modes should be excited when 
a neutron star is formed through gravitational collapse, 
or when two neutron stars merge at the final stage of 
a binary systems evolution, but how much energy goes 
into the pulsation? At the present time we cannot answer 
such questions: To perform the required calculations in 
a relativistic framework is simply beyond our means. 

As a first step towards understanding the issue we have 
studied scattering of axial gravitational wave-packets by 
a compact star. Since axial perturbations are governed 
by a single wave equation with an effective potential 
|p|,p| , this problem is in many respects identical to that 
for black holes [n~3(l . 
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FIG. 1. The response of a uniform density star 
(M/R = 0.2) to a Gaussian pulse of axial gravitational waves. 
The top panel shows the axial perturbation as seen by a dis- 
tant observer, while the lower panel shows the same function 
on a logarithmic scale. 

The result of a typical simulation is shown in Figure 
1. The exponential ringdown at late times (from t « 
100M) corresponds to the first axial w-mode. A power 
spectrum for the data in Figure 1 shows that the first 
three axial modes are excited. A significant part of the 
initial energy clearly goes into quasinormal-mode ringing. 
To study the relation between the fluid modes and the w- 
modes we must turn to the polar problem. Preliminary 
studies pi indicate that i) a polar perturbation typically 
excites both the fluid modes and the w-modes and ii) a 
considerable amount of energy is released through the 
w-modes. 

Our simulations suggest that the w-modes may be rel- 
evant in many dynamical processes involving neutron 
stars. We believe that the modes ought to be excited 



when a stellar core collapses to form a neutron star. Most 
of the initial deformation of spacetime could then be radi- 
ated away in terms of w-modes. For the mode-excitation 
to be considerable the collapse must, of course, be asym- 
metric. To estimate the asymmetries in a realistic su- 
pernova is difficult, but the evidence is that the aver- 
age velocity of radio pulsars is large. Asymmetries in 
the core collapse would explain such high velocities [jl6| . 
Moreover, since it is the asymmetries of the spacetime 
that are relevant for the excitation of the w-modes the 
situation for these modes is not too different from that 
for black holes. It is well-known that both axial and po- 
lar quasinormal modes are generic features of collapses 
that form a black hole flTf] , and that the radiation fol- 
lowing the collision of two black holes is dominated by 
quasinormal-mode ringing flisfl . Hence, it does not seem 
unreasonable to assume that the w-modes will be excited 
when a neutron star is created or when two such stars 
collide. 

One important conclusion that follows from our simple 
simulations is that studies of dynamical processes involv- 
ing neutron stars must be performed using general rela- 
tivity. This may seem obvious, but at the present time a 
considerable effort is invested in calculations using New- 
tonian gravity. The gravitational radiation is typically 
extracted by means of the quadrupole formula. The w- 
modes will never appear in such calculations and thus 
this feature of the full problem will be unaccounted for. 

It is clear that the stellar pulsation modes can play 
a role in many astrophysical scenarios, but will we be 
able to observe them with future gravitational-wave de- 
tectors? In the high-frequency regime where both the 
fluid f-mode and the w-modes reside (1-2 kHz and 8-12 
kHz, respectively) one would not expect too much from 
the new generation of laser-interferometric gravitational- 
wave detectors. These will probably not be sensitive 
enough at such high frequencies (although their perfor- 
mance in this regime may be much improved through 
dual recycling). But neutron stars pulsating in the f- 
mode are ideal sources for the resonant bar detectors that 
are currently operating. The frequency of the f-mode 
also makes it well suited for detection by the recently 
proposed spherical solid-mass detectors, e.g., the Trans- 
verse Icosahedral Gravitational Wave Antenna (TIGA) 
p9| . Meanwhile, the w-modes provide interesting sources 
for the proposed array of smaller bar detectors, which 
should be sensitive in the few kHz regime |2C|j . 

To obtain rough estimates for the typical gravitational- 
wave amplitudes from a pulsating star we use the stan- 
dard relation for the gravitational-wave flux jnj 
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which is valid far away from the star. Combining this 
with i) dE/dt = E/2t where r is the e- folding time of the 
pulsation and E is the available energy ii) the assumption 
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that the signal is monochromatic (with frequency /) and 
iii) the knowledge that the effective amplitude achievable 
after matched filtering scales as the square-root of the 
number of observed cycles: h e g = h^/n = h-s/Jr, we get 
the estimates 



h e ti ~ 3 x 10" 
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for the f-mode, and 

E \ 1/2 /10kHz\ 1/2 /50kpc 
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for the fundamental w-mode. Here we have used typi- 
cal parameters for the pulsation modes, and the distance 
scale used is that to SN1987A. In this volume of space 
one would not expect to see more than one event per ten 
years or so. However, the assumption that the energy 
release in gravitational waves in a supernova is of the 
order of 10~ 6 M c 2 is very conservative 21 1. If a sub- 
stantial fraction of the binding energy of a neutron star 
were released through the pulsation modes we could hope 
to see such events out to the Virgo cluster (and perhaps 
as many as a few per year). 

What can we learn from observations ?. — If they were 
detected, the stellar pulsation modes carry key informa- 
tion about the nature of neutron stars. Suppose that we 
detect a gravitational-wave signal from a compact star, 
what can we hope to learn from it? Again, a detailed 
answer requires much further study. But it is useful to 
speculate on the possibilities. 

Assume that we detect a signal and manage to extract 
both the fluid f-mode and the slowest damped polar (or 
axial) w-mode from it. Then spectral studies suggest 
that: i) the oscillation frequency of the f-mode scales with 
the average density of the star through y/M/R 3 , and ii) 
the damping rate of the w-mode depends linearly on the 
compactness ratio M/R of the star. These properties are 
illustrated in Figures 2a, b. In principle, one should be 
able to infer the average density of the star from an ob- 
servation of the f-mode. Similarly, an observation of a 
w-mode leads to an estimate of the stars compactness. 
Although important, the information available in either 
case would not provide detailed information about the 
stellar parameters. But the situation is markedly differ- 
ent if both modes are observed. Then both the mass and 
the radius of the star can be extracted from the observed 
data. This information puts strong constraints on the 
nuclear equation of state [g2| . 

The idea seems simple enough, but will it be useful 
in practice? Let us give an example for stars with poly- 
tropic equations of state: We have constructed a set of in- 
dependent polytropic stellar models (p — Kp 1+1 / n ) with 



varying polytropic index (n = 0.5; 0.8; 1; 1.2; 1.5). We 
have determined the f-mode and the slowest damped po- 
lar w-mode for each of these models (the results would 
be similar if we use the fundamental axial w-mode) . The 
relevant data is graphed in Figures 2a, b. The theoreti- 
cal spectra are well fitted by the following two relations: 
The oscillation frequency of the f-mode varies with the 
average density of the star as 
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while the damping rate of the w-mode behaves as (r u 
the e-folding time of the pulsation) 



1 , 

— (MHz) 



0.107 - 0.069 



/I0km\ / M 
[~1T 



I AM, 



(5) 



These approximations are shown in Figures 2a, b. The 
crucial question is: How robust are these fits for mani- 
festly independent stellar models? The two relations (Q) 
and (|[) must provide reasonably accurate estimates for 
both M and R for all stars in our dataset in order for the 
idea to be useful. The scheme passes this first simple test 
with flying colours: By inverting (|J) and the mass 
and the radius of each star can be determined at least 
as accurately as can the average density and the stellar 
compactness. The error in M and R is, in fact, typically 
smaller than 10% . 




0.18 0.22 0.26 
compactness 



3 



FIG. 2. (a) The pulsation frequency of the f-mode as a 
function of the average density (in units 10 14 g/cm 3 ) of the 
star , and (b) the e-folding time for the slowest damped polar 
w-mode as a function of the stellar compactness (dimension- 
less). The behaviour is not strongly dependent on the poly- 
tropic model used, and therefore ideal for parameter extrac- 
tion. The polytrope data that is used in the discussion in the 
text is shown as diamonds connected with solid lines. The 
corresponding least-square fits are shown as dashed curves. 
Also indicated (by crosses) are preliminary results for a few 
realistic equations of state. 

This is clear evidence that our idea can be useful in 
a real detection situation, and that it deserves to be in- 
vestigated in more detail. One must: i) incorporate the 
estimated effects of statistical errors and measurement 
ones. It will, for example, be much more difficult to infer 
the w-mode damping rate from a data set than to find the 
f-mode pulsation frequency (remember that the observ- 
ability of a periodic signal buried in noise scales roughly 
as the square- root of the number of observed cycles), ii) 
obtain fits similar to (Ji|) and (||) for more realistic equa- 
tions of state. If such relations prove as robust as the 
ones we have obtained for polytropes then the suggested 
scheme looks truly promising. 

We are presently investigating these questions and will 
return to them in the future. As yet, we have done pre- 
liminary calculations for seven of the equations of state 
that were used by Lindblom and Detweiler || (their mod- 
els A, B, C, E, F, G and I). These results, which have 
been included in Figures 2a, b, are in impressive agree- 
ment with the polytrope data used in our example. As a 
further demonstration of the robustness of the suggested 
scheme we have considered the f-mode data that Cutler, 
Lindblom and Splinter obtained for several realistic equa- 
tions of state |g3| . One can use the stellar parameters M 
and R from Table 2 in |2j| together with our relation (|J) . 
This leads to approximate f-mode frequencies that should 
be compared to the data in Table 3 ofjE3| . Despite being 
derived for polytropes, our equation (14) approximates al- 
most all the £ = 2 frequencies in Table 2 of |23| to well 
within 10 %. 

Concluding remarks. - In this letter we have pre- 
sented suggestive results for problems that can have di- 
rect impact on the future detection of gravitational waves 
from neutron stars. 

First we have shown that the gravitational- wave modes 
of a neutron star will be excited in a dynamical situa- 
tion. The modes should be excited when a neutron star 
is formed in a gravitational collapse, a process when a 
considerable amount of energy is released. But at present 
it is difficult to estimate the amount of energy that goes 
into the pulsation modes. What is absolutely clear (and 
this is a very important conclusion) is that we are asking 
questions that can never be answered within Newtonian 
gravity. The assertions presented here must be tested by 
more detailed, fully general relativistic, simulations. This 



is a challenge for numerical relativity. A relativistic de- 
scription of gravitational collapse to form a neutron star, 
or the merger of two stars, should tell us whether the 
w- modes are of observational relevance or not. Further- 
more, the pulsation properties of a slightly perturbed star 
can be used as a powerful test of the reliability of a nu- 
merical code. 

The second part of our work shows that the stellar pul- 
sation modes can, if observed, provide us with accurate 
information about the stellar parameters. We have il- 
lustrated how a simple scheme can lead to estimates of 
both the mass and the radius of the star. This is a very 
useful suggestion, since it can put strong constraints on 
the nuclear equation of state. 

Given the obvious importance for high-frequency 
gravitational-wave detectors, further studies of the issues 
raised here are urgently required. 
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